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COSMOLOGICAL BLACK HOLE SPIN EVOLUTION BY MERGERS AND ACCRETION 
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ABSTRACT 

Using recent results from numerical relativity simulations of black hole mergers, we revisit previous 
studies of cosmological black hole spin evolution. We show that mergers are very unlikely to yield large 
spins, unless alignment of the spins of the merging holes with the orbital angular momentum is very 
efficient. We analyze the spin evolution in three specific scenarios: (1) spin evolves only through mergers, 
(2) spin evolves through mergers and prolonged accretion episodes, (3) spin evolves through mergers and 
short-lived (chaotic) accretion episodes. We study how different diagnostics can distinguish between these 
evolutionary scenarios, assessing the discriminating power of gravitational- wave measurements and X-ray 
spectroscopy. Gravitational radiation can produce three different types of spin measurements, yielding 
respectively the spins of the two black holes in a binary inspiral prior to merger, the spin of the merger 
remnant (as encoded in the ringdown waves), and the spin of "single" black holes during the extreme 
mass-ratio inspiral (EMRI) of compact objects. The latter spin population is also accessible to iron-line 
measurements. We compute and compare the spin distributions relevant for these different observations. 
If iron-line measurements and gravitational-wave observations of EMRIs only yield dimcnsionlcss spins 
j = J/M 2 > 0.9, then prolonged accretion should be responsible for spin-up, and chaotic accretion sce- 
narios would be very unlikely. If only a fraction of the whole population of low-redshift black holes spins 
rapidly, spin-alignment during binary mergers (rather than prolonged accretion) could be responsible for 
spin-ups. 

Subject headings: cosmology: theory - black holes - gravitational waves - galaxies: evolution 



1. INTRODUCTION 

In general relativity, the properties and dynamics of as- 
trophysical black holes (BHs) are determined only by their 
masses and spins. There is observational evidence that 
BHs with masses in the range ~ 10 6 — 10 9 Mq exist in 
the bulges of nearly all local, massive galaxies, including 
our own. BH mass measurements show an almost-linear 
relation between the BH mass and the mass of the galac- 
tic bulge hosting the B H. Spin measurement s are generally 
more uncertain (|Naravanll2005t lMillerll2007ft . and the cos- 
mological coevolution of BH masses and spins is a signif- 
icant o pen problem in our understanding of quasars and 
AGNs (iGammie et aT1l2004t IShapirdl2005h. 

Massive BH formation scenarios (jVolonteri et al. 2003) 
are important in the planning of the space-based 
gravitational-wave detector LISA. Cosmologi cal BH evo- 
lution affects the rates of detectable event s (|Bertil [2 006: 
Sesana et al.ll2007t) . data analysis strategies (|Arnaud et al.l 
20071) and the science performance of the instrument, in- 
cluding its parameter estimation capabilities. BH spins 
have a strong impact on source modeling and parameter 
estimation, because they induce precession in the orbits of 
inspiralling BH binaries and affect the radiation emitted 
in the merger and ringd own phases. 

I Volonteri et a l. (2005) (henceforth paper I) studied the 
distribution of massive BH spins and its evolution with 
cosmic time under the combined effect of mergers and ac- 
cretion in the context of hierarchical galaxy formation the- 
ories. They found that gas accretion affects the spin evo- 
lution more than mergers, and that prolonged accretion 
efficiently spins holes up to j ' ~ 1. 
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Until recently, the spin of the BH resulting from generic 
BH mergers was est i mated by heuristic argumen t s (see 
e.g. iMerritt k Ekeri (|2002h L iHughes k Blandfordl (|2003h 
(henceforth HB) estimated the merger remnant's spin us- 
ing an extrapolation to comparable-mass binaries of results 
valid for small mass ratios, q — M%/Mi <C 1, showing 
that mergers typically spin BHs down. Their calculations 
become unreliable for "major" mergers with q ~ 1. Re- 
cent breakthroughs in numerical relativity have provided 
us with a quantitative understanding of the spins and re- 
coil velocities resulting from c omparable-mass binary BH 
mergers: see iPretoriusI (|2007f) for a review. Here we up- 
date the analysis of paper I, that used the HB model. For 
the first time, we implement numerical relativity results 
to compute the BH spin resulting from generic mergers in 
hierarchical models of BH formation. We study how the 
efficiency of spin alignment with the orbital angular mo- 
mentum L affects upper and lower limits on the remnant's 
spin by considering three different scenarios: (I) complete 
isotropy, (2) efficient alignment of the spin of the more 
massive hole with L, so that the smaller hole orbits in the 
equatorial plane of the larger, and (3) alignment of both 
spins with L. 

Spin- up is a natural consequence of prolonged disc- mode 
accretion: any hole that has doubled its mass by capturing 
material with consta nt angular momentum axis w ould end 
up spinning rapidly (|Bardeenllf970l : lThorndlf974t ). How- 
ever, when the angular momentum of the accretion disc is 
misaligned with respect to the direction of the BH spin, ac- 
cretion of counter-rotating material can cause a spin-down 
of the hole. The orbits in the inner accretion disc co-rotate 

, CA 91109, USA 
USA 



2 



Berti & Volonteri 



or counter-rotate depending on the ratio between the an- 
gular momentum of the disc and of the hole: if the cosine 
of the inclination angle cos</> < — Jd/(2j), where jd is the 
total angular momen tum of the accre tion disc, then the 
disc counter-rotates (|King et al.ll2005l ). Sustained accre- 
tion from a counter-rotating disc spins down a maximally- 
rotating hole to j = after the hole has increased in mass 
by a factor ~ a/3/2 (|Bardeenl ll970). A complete overflip 
eventually occurs, and then accretion of co-rotating mate- 
rial acts to spin up the BH again: a 180° flip of the spin 
of an extreme-Kerr hole will occur after the BH triples in 
mass. Paper I argued that the lifetime of quasars is long 
enough that the innermost regions of accretion discs align 
with BH spins (possibly through spin flips), and hence all 
AGN BHs should have large spins. 

The picture of s pin evolutio n via p rol onged accretion 
was q uestioned by IKing et all (|2005l ). iKing fe Pringld 
(2006i) suggest that accretion always proceeds via very 
small (and short) episodes, caused by fragmentation of 
the accretion disc where it becomes self-gravitating. The 
hole would then accrete each time an amount of mass 
corresponding to a tiny fraction of its own mass, and all 
episodes would have uncorrelated directions ("chaotic ac- 
cretion"). Since counter-rotating material spins BHs down 
more efficiently than co-rotating material spins them up, 
this scenario implies that BH spins are very small: ac- 
cretion of randomly oriented droplets of gas would rapidly 
spin down any BH that had its spin increased by a "major" 
merger. The only rapidly spinning holes would then be 
those which have recently experienced a merger, and have 
not accreted any matter afterwards. IKing et all |2008) 
refined their chaotic accretion scenario by statistical con- 
siderations, concluding that accretion very rapidly adjusts 
the hole's spin parameter to average values j ~ 0.1 — 0.3 
from any initial conditions, but with significant fluctua- 
tions Aj ~ ±0.2 about these values. In this paper we 
compare the prolonged-accretion scenario and the chaotic 
accretion scenario, showing that rapidly spinning holes 
should be extremely rare if chaotic accretion is the norm. 

In Section [2] we review numerical simulations of BH 
mergers and discuss their predictions for the final spin. In 
Section [3] we present our main results on cosmological spin 
evolution in different merger scenarios, with and without 
accretion. 

2. FINAL SPIN FROM BINARY MERGER SIMULATIONS 

In the last two years there has been enormous progress 
in the numerical simulation of BH binaries. We now know 
that equal-mass non-spinning binaries produce a rotat- 
ing (Kerr) BH with final dimensionless spin parameter 
|j fin | ~ 0.69. For unequal-mass non-spinning mergers with 
mass ratio q, the final spin is well fitted by the sum of 
two terms: |j fin | ~ 2 v / 3q/(l + q) 2 - 2.029g 2 /(i + q) 4 . The 
first term is an extrapolation to comparable masses of the 
orbital angular momentum of a particle at the innermost 
stable circular orbit of a non-rotating BH, and the second 
term account s for the angular m omentum radiated in the 
final plunge (jBerti et al.ll2007al) . 

A significant sample of spinning binary BH merger 
simulations is now available. Different groups showed 
that when the initial spins are large and aligned with 
the orbital angular momentum L the binary "hangs 



up" , radiating more energy and angular momentum 



(ICam panelli et al. _ 200; 



Bak er et all l2007t iPollnev et all 
2007: iMarronetti et al] I2007D . Even accounting for the 
additional angular momentum due to orbital eccentric- 
ity, no violation of cosmic censorship should occur: a 
binary BH merger should always result in the forma- 
tion of a Kerr BH (jSperhake et al.ll2007f ). Simulations 
have also been carried out for some "generic" orientations 
of the initial spins, including configur ations leading to 
spin- fl ips or to Schwarzschild r emnants (ICampanelli et"ail 
[20071: iHerrmann et alj l2007t iTichv fe Marronettil l2007h 
iBerti et al.ll2007bD . 

Unfortunately, covering the whole parameter space by 
numerical simulations is computationally costly and not 
practical. Semi-analytical models for the Kerr parameter 
of the final hole are particularly useful for astrophysical ap- 
plications. One such model , based on point-particl e analo- 
gies, has been proposed bv lBuonanno et aTT(|2007l ). Their 
model can reproduce the final spin computed by numeri- 
cal simulations within a few percent, the disagreement be- 
ing larger when both spins are large and anti-aligned with 
L. Here we use semi-a nalytical fitting formul as of numer- 
ical results derived bv lRezzolla et al.l (|2007a .b ,c). Unfor- 
tunately, due to the inaccuracy of large-spin simulations 
and to errors in the fitting parameters, the final spin |j fin | 
predicted by the fits can be slightly larger than the Kerr 
bound when the initial Kerr parameters of the holes |jj 
(i = 1, 2) are large and aligned with L. To get rid of 
this undesired feature we simply introduce a cutoff on the 
final spin at |j fin | = 1. An elegant and potentially very 
accurate recipe to compute the final spins has been pro- 
posed by iBovle et all (|2007l) using symmetry arguments. 
Some parameters in their "spin expansion" are presently 
undetermined, but in principle they can be fixed by a rea- 
sonably small number of dedicated numerical simulations. 

The spin expansion model should be able to provide 
accurate and general predictions for the final spin in 
the nea r future, but the fitting formula of iRezzolla et alj 
(2007a) is accurate enough for our preliminary ex- 
ploration. Their formula can be used to compute 
|jfinl(<7> I J 1 1 ,\fo\, cosa, cos/3, cos 7), where cosa = jj -j 2 , 
cos/3 = j± • L and cos 7 = j 2 ■ L. 

As a starting point, for comparison with HB, we com- 
puted the final spin resulting from a merger where the 
smaller BH is non-spinning (|j 2 | = 0). Since |j 2 | = the 
angles cosa and cos 7 are irrelevant, and only cos/3, the 
orbital inclination of the smaller hole, matters. Contour 
plots of the final spin |j fin | in the (cos/3 , jj-J) plane, as ob- 
tained from the fitting formula, are shown in our Figure [TJ 
that should be compared with Figure 1 in HB (notice that 
our angle (3 is denoted by 1 in their work). For q = 1/10 
our calculations are in nice agreement with HB: this is a 
useful consistency check of our approach in the small mass- 
ratio limit. However, observational and theoretical argu- 
ments suggest that the coalesc ence of comparable-mas s 
BHs should be rather common (|Bogdanovic et aI1l2007l ). 
When q ~ 1 and |j 2 | = 0, we find that the range of vari- 
ability of the final spin is sensibly smaller than predicted 
by HB: 0.3 < |j fi J < 0.9 for q = 1/2, and 0.5 < |j fin | < 0.8 
for q = 1. 

In Figure [2] we drop the assumption |j 2 | = and we 
show contour plots of |j fin | in the , |j 2 |) plane, for se- 
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lected values of q. To deal with the angular dependence 
we consider three different merger scenarios: 

(1) Isotropic mergers: for each pair of initial spin mag- 
nitudes (|j]J , |j 2 |) we average over (cos a , cos (3 , cos 7), as- 
suming isotropy on all three angles. This situation should 
be common in "dry" mergers, i.e. when holes do not ac- 
crete during merger, evolving solely via stellar dynami- 
cal processes. In vacuum, post-Newtonian equations pre- 
dict that spin-orbit resonances will produce alignment 
of the spins in a very small region of par ameter space 
(|Schnittmanl 12004 IBogdanovic et all |2007), so isotropy 
should be a good assumption in the absence of accretion 
discs or gas. For all values of q, the minimum final spin 
results from non-spinning mergers with Ij-J = |j 2 | = 0. 
Naively we could expect that the average over all angles 
(cos a , cos (3 , cos 7) should produce a final spin that is very 
close to the spin produced by non-spinning mergers, at 
least for equal masses. As we consider larger mass ratios 
the larger BH should play a more important role, so (on 
average) the final spin should be slightly larger than the 
value predicted by non-spinning mergers. This expecta- 
tion is confirmed by the top row of Figure [21 In general, 
our calculations are in agreement with the HB conclusion 
that isotropic mergers tend to "spin-down" a fast-spinning 
hole. 

(2) (Anti) aligned mergers: cos a = ±cos/3 = ±C0S7 = 
1, so the BH spins are always parallel (aligned or an- 
tiali gned) with L. In the " wet merger" scenario proposed 
by IBogdanovic et"al] (|2007l ). alignment (upper signs in ±) 
should be more likely than antialignment (lower signs), 
because the cumulative angular momentum of the ac- 
cretion disc is much larger than the angular momentum 
of th e BHs (but see King et all 120051 : lLodato k Pringld 
2006). If alignment is indeed the norm, we should focus 
on the top right quadrant of the contour plots in the cen- 
tral row. Then the minimum spin is still obtained when 
liil — I j 2 1 = 0, as in the isotropic case. The main difference 
is that now the maximum spin can be quite close to one, 
with a lower limit |j fln | ~ 0.96 for equal-mass mergers. 

(3) Equatorial mergers: this case is intermediate between 
cases (1) and (2). Here the smaller BH inspirals in the 
equatorial plane of the larger (cos (3 = 1), but for each pair 
(Ijili U2I) we average assuming an isotropic distribution in 
(cos a, cos 7). Equatorial mergers could occur if Newto- 
nian dynamical friction in a flattened system brings the 
smaller hole into the plane of a central (gaseous or stellar) 
disc before it gets into the relativistic regime. All mergers 
with q > 1/2 produce an average final spin |j fin | < 0.92, 
and |j fin | < 0.85 for q = 1. As in case (1), to get large 
Kerr parameters |j fin | > 0.9 in the absence of prolonged 
accretion we must start with large spins and avoid major 
mergers. 

For each value of q, it is interesting to maximize and 
minimize the final spin resulting from a merger in the 
(Ijilib'al) pl ane - The resulting extrema in the three dif- 
ferent scenarios are plotted as functions of q in Figure [3] 
For each value of q we average over the individual spin 
magnitudes C C I J 1 1 ^ I J 2 1 ) - Both (j-J and |j 2 | are assumed to 
be uniformly distributed between and 1 . In the isotropic 
case, the allowed average values of the final spin range 
between the solid and dashed black lines. The most in- 
teresting feature is the appearance of a funnel at q > 0.1: 



major mergers tend to produce BHs with average spins 
very close to the value |j fin | ~ 0.69 resulting from equal- 
mass, non-spinning mergers. For smaller q the larger BH 
dominates the dynamics, and the final spin can be sub- 
stantially larger or smaller than this value. Moving from 
isotropic to equatorial and aligned mergers can produce 
slightly larger maximum spins (as indicated by the solid 
lines of decreasing thickness) , but in all three scenarios the 
most likely spin resulting from "major" mergers is very 
close to |j fin | ~ 0.69. We now explore some consequences 
of these results for spin evolution in cosmological BH for- 
mation scenarios. 

3. COSMOLOGICAL SPIN EVOLUTION 

In Figure 3] we show histograms of the spin distribution 
of merging BHs in different redshift ranges. For simplicity, 
BH seeds at high redshift are assumed to be non-spinning. 
These plots are useful to isolate the contribution of merg- 
ers to the spin evolution of the whole BH population, and 
they are of direct interest for LISA obs ervations of mas- 
sive BH inspiral, merger and ringdown teerti et al]|2005L 
l2006t lLang k Hughesj|2006t iKocsis et alj|2007j ). The three 
columns in each of the three plots correspond to the merg- 
ing scenarios described above: isotropy, efficient alignment 
of the spins with L, and equatorial mergers. 

In the left plot, spin evolution is due to mergers only 
and accretion is ignored. Consider first isotropic mergers 
(left column). At high z most BH seeds have compara- 
ble masses, and since mergers with q > 0.1 produce rem- 
nants with |j fin | ~ 0.7 (Figure[3]) the spin distribution post- 
merger peaks around this value. Later on, small-q mergers 
become more common, and on average (as shown in HB) 
they tend to spin down the hole, so the spin distribution 
flattens out. For aligned mergers (central column) the evo- 
lution is different. At high z, mergers of comparable-mass, 
non-spinning holes again produce a peak at |j fin | ~ 0.7. 
However, if alignment is efficient mergers of spinning holes 
have a tendency to spin up the remnant, and eventually 
most BHs at small z are rapidly spinning. In the equatorial 
merger scenario (right column) the less massive hole has 
marginal impact on the dynamics, and the overall spin dis- 
tribution is qualitatively very similar to the case of align- 
ment. 

In the central plot spin evolution is due to both merg- 
ers and accretion, where accretion is modeled as in paper 
I: prograde or retrograde orbits are equa lly probable, and 
spins evolve according to iBardeenl (|1970f ) . If enough mass 
is available to the hole, and the disc was initially counter- 
rotating with respect to the hole, an overflip can eventually 
occur. In this scenario, accretion-induced spin-up is very 
efficient. For aligned (equatorial) mergers, > 90% of merg- 
ing BHs at z < 5 [z < 2) have j > 0.9. If mergers occur 
isotropically, on average they tend to spin down the holes, 
partially counteracting the accretion-induced spin-up. As 
a result, in this case the spin distribution post-merger has 
a long tail extending all the way down to j ~ 0. 

In the right plot we consider the chaotic accretion sce- 
nario. Now accretion happens in short-lived episodes, 
where the BH increases its mass by 0.1%, and prograde 
or retrograde orbits are equally probable. Mergers per 
se are unlikely to produce fast-spinning remnants, and 
chaotic accretion is very efficient in spinning BHs down. 
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Therefore, at all but the highest redshifts the spin dis- 
tribution post-merger is roughly uniform in the range 
j G (0, 0.7 — 0.8), and the chance of mergers involving 
high-spinning BHs is very low. LISA measurements of 
spins j > 0.9 would strongly favor prolonged accretion 
and indicate that alignment is efficient in merger events. 
Conversely, measuring values of j < 0.9 at z < 5 would 
indicate that accretion is chaotic (or negligible), or that 
alignment is not efficient. 

For BH sp in measurements based on X-ray spectroscopy 
(Miller 2007) and for the low-redshift observations of gravi- 
tational waves from compact objects falling into a massive 
BH (EMRIs) , that could allow ve ry accurate spin mea- 
surements ()Barack fe Cutler] l2004f ) , we are interested in 
the spin distribution of the whole BH population. This 
distribution is shown in Figure [5j 

When we consider the spin evolution under mergers only 
(left plot) the overall spin distribution is bimodal. One 
peak is around j ~ because we assume that seed BHs 
are non-spinning, and some BHs never experience merg- 
ers. A second peak is located around j ~ 0.7 (for isotropic 
mergers) and it extends all the way up to j ~ 1 for 
aligned/equatorial mergers. 

If accretion influences the evolution of BH spins, it dom- 
inates over mergers in determining the spin evolution of 
the whole BH population, confirming the results of Paper 
I. We stress that our accretion model is highly oversimpli- 
fied, and does not take into account the different evolution 
of spins in different galactic environments: here we sim- 
ply assume that accretion is triggered always and only by 
galaxy mergers. This is indeed a simplisti c assumption, 
especi ally for low-redshift faint AGNs: see I Volonteri et all 
(|2007| ) for a detailed analysis of the connection between 
spin evolution and galaxy morphology. 

Whatever the dominant merger scenario, under pro- 
longed accretion (central plot) spin- up is very efficient, and 
a large fraction of all BHs has spin j > 0.9. On average 
isotropic mergers tend to spin BHs down: if isotropy is the 
norm, ~ 30% of the whole BH population has j < 0.9, with 
a roughly uniform distribution. For aligned and equatorial 
mergers spin-down is less efficient, and only ~ 10% of all 
BHs have j < 0.9. 

The chaotic accretion scenario (right plot) predicts a 
completely different distribution. In this case most BHs 
(~ 50 — 80% of the total) have j < 0.1. When the average 
spins are so low mergers act to spin BHs up, and the spin- 
up is more efficient in the aligned and equatorial cases. 
This produces a roughly uniform spin distribution in the 
range j £ (0.1, 0.7). Interestingly, the tail of the distribu- 
tion never extends above j ~ 0.8, even for small redshifts. 
The measurement of lar ge BH spins, such as th e value 
of j ~ 0.99 claimed by iBrenneman fe Reynolds! 12 006) 
for M CG-06-30-15 (see also lElvis et al.ll2002t IWang et al l 
2006) would presumably indicate that chaotic accretion is 
not the norm. In any case, BH spin measurements should 
easily distinguish between chaotic and prolonged accre- 
tion. Both for merging BHs and for the whole BH pop- 
ulation, under chaotic accretion maximal Kerr BHs are 
extremely unlikely, and the distribution peaks at j < 0.7 
for all z < 10. 

Suppose that LISA measurements of the spin resulting 
from a binary merger yield j > 0.9. Then it is hard to 



tell if spin is accretion-independent, but alignment is ef- 
ficient (left plot in Figure [4j or if instead prolonged ac- 
cretion is responsible for spinning BHs up (central plot). 
Luckily, this degeneracy is broken when we consider the 
whole population. If accretion is negligible and mergers are 
not isotropic, only < 50% of all BHs have spins j > 0.9. 
If instead spin-up is due to accretion, the population is 
largely dominated by spins j ~ 1 (compare the left and 
central plots in Figure©. EMRIs (jBarack fc Cutleij|2004h 
and X-ray observations of the Fc Ka line b y planned mis- 
sions such as Constellation-X (Miller 2007) will probe the 
whole population at low redshifts (z < 2), while observa- 
tions of the inspir al (|Lang fe Hughes! l2006f ) and ringdown 
(jBerti et al.ll2006f) of massive BHs can provide a census of 
BH spins over a wider range of epochs. A combination of 
these techniques (and possibly others) has the potential to 
probe unequivocally the physical processes involved in the 
evolution of BH spins. 

4. CONCLUSIONS 

Understanding how fast BHs spin is as important as un- 
derstanding how they grow in mass. The spin and mass 
evolution of BHs are intimately linked. The expected spin 
of a hole depends on whether it gained most of its mass 
via mergers or accretion. Conversely, the spin influences 
how efficiently BHs accrete mass, determining the mass- 
to-energy conversion efficiency in radiatively efficient ac- 
cretion phases. The spin also determines how much energy 
can be extracted from aBH. Assuming that relativistic 
jets are powered by rotating BHs through the Blandford- 
Znajek mechanism, the so-called "spin paradigm" asserts 
that powerful rela tivistic jets are p roduced in AGN with 
fast rotating BHs (|Blandfordll99of ). 

Studying the role of BH mergers in the evolution of 
BH spins has been daunting until recently. Due to re- 
markable advances in numerical relativity, we finally have 
a quantitative understanding of the spins result ing from 
comp arable-mass binary BH mergers (see e.g. iPretoriui 
12007t h Here we use numerical relativity results to investi- 
gate plausible scenarios of spin evolution through mergers 
and accretion in a cosmological context. The co-evolution 
of BH masses and spins is studied self-consistently: merg- 
ers and accretion determine BH spins, which in turn lend 
to the calculation of the radiative efficiency during accre- 
tion episodes. 

We focus on three scenarios for the mass and spin co- 
evolution: (1) spins evolve only through mergers, (2) spins 
evolve through mergers and prolonged accretion episodes, 
(3) spins evolve through mergers and short-lived (chaotic) 
accretion episodes. If BHs accreted most of their mass 
through prolonged disc-mode accretion, by adding mate- 
rial with constant angular momentum axis, they would end 
up spinning rapidly. If instead BHs built-up their mass via 
short-lived episodes with uncorrelated angular momentum 
axis, the typical spin of BHs would be very low. 

We further consider how the dynamics of BH mergers in- 
fluences the final spin. The mutual directions of the spins 
of the holes in a binary with respect to the orbital plane 
strongly affect the final spin of the remnant. If the spins 
align efficiently with the orbital angular momentum, the 
spin of the merger remnant is larger than in the case of 
random orientations (see Figures [H and [3]) . We consider 



Cosmological black hole spin evolution 



5 



three physically motivated cases: (1) random orientations 
(isotropy), which is expected if BH binaries are in gas- 
poor host galaxies, or if the orbital plane is not aligned 
with the disc plane; (2) alignme nt between spins and or - 
bital angular momentum, which [Bogdanovic et alj (|2007t ) 
suggest should be the norm for binaries forming during 
gas-rich galaxy mergers; (3) equatorial orbits, where the 
smaller BH inspirals in the equatorial plane of the larger. 
Equatorial mergers could occur if dynamical friction brings 
the smaller hole into the plane of a central disc before it 
gets into the relativistic regime. 

Except in the case of aligned mergers, we find that a 
sequence of BH mergers can lead to large spins |j fin | > 0.9 
only if BHs start already with large spins and they do 
not experience many major mergers. This is illustrated 
in Figure [3J where we show that (on average) isotropic 
mergers with q > 0.05 lead to |j fin | < 0.9. Therefore, the 
common assumption that mergers betw een BHs of similar 
mass always lead to large spins (e.g., IWilson k, Colbert! 
Il995h needs to be revised. In the isotropic case, numer- 
ical relativity results imply that major mergers tend to 
produce BHs with average spins very close to the value 
|j fin | ~ 0.69 resulting from equal-mass, non-spinning merg- 
ers (see again Figure [3J . 

Our models can be used to test the discriminating 
power of direct observational techniques: gravitational- 
wave measurements and X-ray spectroscopy. LISA mea- 
surements can provide information on the spins of the two 



BHs in a binary prior to merger (inspiral) and on the spin 
of the merger remnant (ringdown). Such measurements 
can unequivocally inform us on the typical spin of the BH 
population (Figured]). However, if LISA measurements of 
the spin resulting from a binary merger yield j > 0.9 a de- 
generacy remains: a distribution skewed towards large val- 
ues can occur if spin is accretion-independent, but align- 
ment is efficient (left plot in Figured]), or if instead pro- 
longed accretion is responsible for spin up (central plot). 

This degeneracy is broken by coupling spin measure- 
ments of binaries to either X-ray spectroscopy or EMRIs, 
as these techniques sample the whole population. If iron- 
line measurements and gravitational-wave observations of 
extreme mass-ratio inspirals (EMRIs) only yield dimen- 
sionless spins j = J/M 2 > 0.9, then prolonged accretion 
should be responsible for spin-up, and chaotic accretion 
scenarios would be very unlikely. If instead only a fraction 
of the whole population of low-redshift BHs spins rapidly, 
spin-alignment during binary mergers (rather than pro- 
longed accretion) could be responsible for spin-ups. 
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Fig. 1. — Update of Figure 1 in[Hjighes_& Blandford (2003) . Here we assume |j 2 | = (so that cos a and cos 7 are irrelevant) and we show 
the influence of the orbital inclination cos on the final spin for different values of 
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Fig. 2. — Contour plots of the final spin for different mass ratios (left to right: q = 1, q = 1/2 and q = 1/10, respectively) in the (|jj| , |j 2 |) 
plane. Recall that M2 < Mi. We consider three different scenarios. In the top row we average over the sky assuming isotropy on all three 
angles. In the middle row we assume that alignment is efficient, as proposed by Bogdanovic ct al. (2007). By convention, in this case ji < 
means that the spin of BH i is antialigned (rather than aligned) with the orbital angular momentum. In the bottom row we set cos j3 = 1 (so 
the spin of more massive BH is aligned with the orbital angular momentum, and the smaller hole orbits in the equatorial plane of the larger) 
and we assume isotropy in (cos o , cos 7) . 
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spin for equatorial inspirals. Solid blue (thin) line: maximum spin in the aligned case. Dashed black (thick) line: minimum average spin in 
the isotropic case. Since the minimum is attained when [j-J = |j 2 | = 0, this line is also the minimum average spin attainable by equatorial 
inspirals, or the minimum spin attainable in the "aligned" case when we rule out the possibility of antialignment. Dashed blue (thin) line: if 
we do allow for antialignment, the minimum spin can become negative (we have a spin flip) when q < 0.3 or so. The horizontal (green) line 
corresponds to a final spin |j fin | = 0.9. When q is close to one, such large spins are only achievable if alignment is efficient. 
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Fig. 4. — Spin evolution of merging BHs due to: mergers only (left); mergers and prolonged accretion (center); mergers and chaotic accretion 
(right). In each plot we consider three representative merger scenarios (see text) and we show histograms of the spin distribution for different 
ranges of variability of the redshift z. Hatched (red) histogram: spins of the binary members before merger. Thick (black) histogram: spins 
of BHs after merger. 
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Fig. 5. — Spin evolution of all BHs due to mergers only (left), mergers plus standard accretion (center) and mergers plus chaotic accretion 
(right). 



